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Perturbative estimates suggest that extended topological defects such as cosmic strings emit 
few particles, but numerical simulations of the fields from which they are constructed suggest the 
opposite. In this paper we study the decay of the two-dimensional prototype of strings, domain walls 
in a simple scalar theory, solving the underlying quantum field theory in the Hartree approximation. 
We conclude that including the quantum effects makes the picture clear: the defects do not directly 
transform into particles, but there is a non-perturbative channel to microscopic classical structures 
in the form of propagating waves and persistent localised oscillations, which operates over a huge 
separation of scales. When quantum effects are included, the microscopic classical structures can 
decay into particles. 
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I. INTRODUCTION 

Our current fundamental theory of the universe sug- 
gests higher order symmetries which are successively bro- 
ken as the universe cools. This naturally leads to the 
formation of topological defect networks 1|. Such ex- 
tended objects could play a role in the early formation of 
structure in the universe 0,0, El- Indeed, recent calcula- 
tions [H show that models with cosmic strings and other 
topological defects fit the Cosmic Microwave Background 
data better than the standard power-law ACDM model. 

Analytical estimates suggest that particle production 
from decaying strings (or, in two dimensions: domain 
walls) is suppressed by the separation of cosmological (£) 
and microscopic (M) scales. The former is the curvature 
scale of the string, assumed to be of order the Hubble 
radius, and the latter is determined by the underlying 
particle physics, typically on the GUT scale. Particle 
production would mean transforming the energy in the 
deep infrared into ultraviolet excitations, over something 
like 58 orders of magnitude in momentum scale. If the 
cosmic string network is prohibited from losing energy, 
its dynamics follows the Nambu-Goto action [3, 6]. Per- 
turbation theory suggest that the radiative energy loss is 
~ Mjl [7|. With more insight into the non-linear domain 
wall dynamics one discovers another channel through 
cusp annihilation Q, at a rate of ~ M 5 / 3 ^" 1 / 3 . Once 
one takes into account the gravitational channel a simple 
estimate gives ~ 10 2 M 4 / M 2 , where M p is the Plack mass 
[l| , which therefore seems to be dominant for sufficiently 
large I. 

Numerical analysis of the Nambu-Goto action confirms 
[ToL [Til . OjJ the analytical scaling assumption [l3j] that 
implies a string density M/t 1 (in d = 2, 3 space di- 
mensions), at least for strings as long as the horizon 
size. An exact Minkowski space simulation of a string 



network reveals that fragmentation to loops is the domi- 
nant decay process 0, EH] , which proceeds down to the 
smallest scale on the network, which may even be the 
microscopic scale of the string width [l6| . Simulations 
in an expanding universe [Til . [T2| broadly confirm this 
picture, although indicate that the relevant small scale is 
the initial correlation length [TJ d, Qj| 

However, for strings which are topological defects, we 
can check this picture by solving the underlying field the- 
ory in the classical approximation. This means in prac- 
tice to integrate the non-linear wave equations on a spa- 
tial lattice and to average over an initial ensemble. This 
approach offers a full insight into all non-perturbative 
phenomena, although it may be difficult to justify the 
omission of quantum effects at the microscopic scale. 
Nevertheless, this method has been successfully used else- 
where, e.g. to explore the dynamics of symmetry break- 
ing [2(| , or of non-thermal phase transitions [21( in the 
post-inflationary Universe. 

The scaling behaviour suggested by Nambu-Goto sim- 
ulations is manifest in the classical field dynamics, even 
though only the microscopic scale appears in the equation 
of motion. It has been demonstrated in the context of 
gauge strings in the Abelian Higgs model [22|, HH global 
strings [23j,|24(, non- Abeli an g lobal strings with junctions 
25 1, and semilocal strings [261. 27\, as well as domain walls 
28, [2^|, including models with junctions [30]. The scal- 
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ing is present in Minkowski as well as expanding space 
time, and in two as well as three dimensions. It seems 
to be a universal feature of classical field theories with 
extended structures. 

For strings, a major difference between the numerical 
solutions of the classical field dynamics and Nambu-Goto 
simulations is that defects decay into classical radiation 
[22| . at a much faster rate than anticipated from per- 
turbation theory and cusp annihilation. One typically 
observes a length of string I in a volume £ d , and hence 
that the string length density is L ~ l x ~ d . Given a mass 
per unit length of /i ~ M 2 , the energy density in string is 
M 2 l 1 ~ d . Since scaling implies l^t, the energy loss rate 
per unit length M 2 /£. Hence, for a loop of size I, the 
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average energy loss rate is M 2 , which is, in fact, greater 
than the gravitational estimate. 

We discover then that the classical scaling implies a 
strong radiative decay. This is very puzzling in view of 
the scale separation between the I and M which grows as 
the simulation proceeds. However, apart from confirming 
the scaling over more than three orders of magnitude, it 
is not our purpose to address this important question. 
Instead we note that up to now it has been unclear if 
the classical approach is valid here, where dynamics is 
driven by an interplay between macroscopic and micro- 
scopic scales. A check for quantum corrections is crucial. 

An alternative to the classical approach is the 
two-particle-irreducible (2PI) effective action technique, 
which is based on a selective resummation of pertur- 
bative diagrams [3l| . Preheatin g d ynamics with non- 
perturbative particle production [32| and particle ther- 
malisation by scattering [33j are within the range of its 
applicability. The so far used homogeneous version of this 
elaborate technique is, however, incapable of addressing 
the question of defect formation (34| • 

We can combine these techniques, using the classical 
approach to form defects and then studying their evolu- 
tion in the 2PI framework. If we keep the next-to- leading 
order diagrams in the 2PI effective action, we will gain 
insight into the scattering and thermalisation of the pro- 
duced particles. The inhomogeneous variant of the 2PI 
equations, however, is technically hardly feasible. Keep- 
ing the lowest order 2PI diagram yields an approximation 
scheme, that is equivalent to the well-known Hartree ap- 
proximation [35|, [36| . While scattering between the pro- 
duced particles is not included here, even the homoge- 
neous version of this scheme could account for the non- 
perturbatively rapid particle production in the early Uni- 
verse [I?], HH • The extension of the equations to inhomo- 
geneous backgrounds was historically motivated by the 
hope that the background field could mediate interaction 
between the freely streaming particles. Although numer- 
ics have shown that the opposite was true [H, E3, EH , 
this method can be still used for finding the leading quan- 
tum corrections to the evolution of classical structures, 
as has been suggested by a one-dimensional analysis of 
moving kings [42j. 

In this paper we analyse the classical solution of the 
A3> 4 theory in two space dimensions corrected by the 
Hartree approximation. In the broken phase this toy 
model features domain walls, which resemble strings in 
this low dimensional setting. We check if there is a sig- 
nificant alteration to the kink dynamics by the inclusion 
of this type of quantum correction. 

First we recall the results from classical simulations 
and demonstrate the scaling behaviour also found in 
Ref. [28[ . Then in section IIIII we introduce the Hartree 
approximation of the considered model. Next, in section 
IIVI we numerically compute the domain wall evolution 
both in the classical and in the Hartree approximated 
framework. We discuss possible interpretations of the 
results in section [V] and finally conclude in section IVT1 



II. CLASSICAL DECAY OF DOMAIN WALLS 

A. Model details 

The Lagrangian density of our scalar theory is as sim- 
ple as 

£ = \ m 2 \m*# - y/ (1) 

The theory has a Z(2) symmetry, this breaks sponta- 
neously if the thermal mass turns negative. By the choice 
of the bare mass parameter m we make sure that the sys- 
tem is deeply in the broken phase at zero temperature. 
The used quartic potential is motivated by the simple 
form of the classical kink solution: 

<J>(x, y) = v tanh (Mx) (2) 

with 

M 2 = -to 2 / 2 and v = V-6m 2 /A. (3) 

The tension of the wall is inversely proportional to the 
coupling: a = 4|to| 2 /3A. In the classical limit the actual 
magnitude of the coupling is irrelevant as it can be scaled 
out. In the numerics we used A = 6M. The only other 
parameter, the mass sets the scale for the evolution, we 
use the inverse wall width M to render all variables di- 
mensionless, this numerically means M = 1. 

We discretise the model on a spatial lattice. We solve 
a cut-off theory with a pre-set lattice spacing a. Since 
much of the physics of our interest is in the infrared, a 
plays little role. Based on earlier numerical experience we 
can use lattices as coarse as aM = 0.5. We repeated the 
presented numerical analysis on a coarser lattice (aM = 
0.7) and found no significant difference. The lattice size 
L, however, matters. In order to avoid the interaction of 
a pair of signals originating from the same site we stop the 
simulation at t = L/2. This assumes that at t = there 
is no correlation between any of the sites. Our initial 
condition will approximately satisfy this condition. 

The initial condition can introduce other scales. We 
start the dynamics from a low energy density random 
configuration with a rich domain wall structure. Our 
main interest is how these walls evaporate under the re- 
alistic assumption that the scale in the initial condition 
separates from the microscopic scale M. 

We designed the following numerical experiment. We 
start from a white noise configuration at t = 0, deep in 
the symmetric phase. We also checked the invariance of 
our results under replacing the initial noise by a (classi- 
cal) thermal equilibrium of the same energy density. We 
then apply a cooling by adding a friction term to the ki- 
netic term in the equation of motion: d 2 tfi — > d 2 4>+2^do4>. 
This evolution is non-physical, and we switch off at a con- 
veniently chosen time when the particle content is neg- 
ligible and the domain wall density reached a desirable 
value [53| . We starting the numerical observations only 
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after a short period of relaxation after the non-physical 
dynamics has been switched off. This corresponds to the 
pre-thermalisation time scale [43[ • 

The field configuration at this instant is the initial con- 
dition of the dynamics of interest. We could set the origin 
of time to this instant, but we choose not to. The t = 
point marks the onset of cooling, because at that point 
the correlation length is known to be of the microscopic 
scale. 



B. Scaling solutions 

The solution of classical dynamics is a straightforward 
computational task and being restricted to two spatial 
dimensions our resources allows for larger lattices (L > 
1000), too. 

As a first observation, we can estimate the domain wall 
length following the techniques detailed in Ref. [28j], and 
reproduced their scaling solutions. In Fig. [T] we show the 
inverse domain wall density as a function of time. If the 
length scale (£) of the domain wall network decouples 
from the microscopic scale, one may expect from dimen- 
sional reasons: i~ x ~ L ~ i -1 . The domain wall density 
L we define as the total length of domain walls on the 
lattice divided by the volume. A link on a lattice is con- 
sidered as a part of a domain wall if the sites at its both 
ends have field amplitudes of opposite sign. 
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FIG. 1: Inverse domain wall density as a function of time. 
There is a perfect linear correspondence over at more than 
two orders of magnitude {t = 10 ... 2000), from the time of 
ending the damped evolution at t = 10 to the end of the 
simulation, when the walls' mean curvature radius is over 10 3 
times larger than their width. 

As always in classical field theory we always display an 
average over an ensemble of runs. (This is the thermal or 
white noise ensemble at t = 0. At positive times there is 
no randomness in the dynamics.) The averaged domain 
wall densities start deviating near t f» L, slightly later 
than expected. 

This scaling is a manifestation of a more generic feature 
in classical field theories, as it has been found in flat or 



curved space-time, and in two or three dimensions [28|. 

So that we gain more insight into the observed scaling 
we show a pair of lattice configurations in Fig. [5] The 
similarly looking snapshots were taken at times 50 and 
100, respectively, but the earlier configuration we halved 
in linear size and scaled up accordingly. The time evolu- 
tion appears to be equivalent to zooming. 




FIG. 2: Two snapshots of the lattice field configurations 
(L — 2000). Blue and orange regions show the domains of 
the degenerate vacua. Fo the left, we show the configuration 
at t = 50, we cropped the region (0, L/2) x (0, L/2) and scaled 
up by a factor of two. To the right we show the uncropped 
lattice at t = 100. There is no qualitative difference between 
the snapshots. 

This feature can be made more formal in terms of the 
correlation function. We define our C(r, t) correlation 
function as 

C{r, t) = J dxdydz {cj)(x, z, t)<f>(y, z + r, t)) , (4) 

where (•) denotes an ensemble averaging. If the more 
conventionally defined correlation function [Hij G(\x — 
y\,t) = (cj)(x,t)(j)(y,t)) scales as G(r,i) = t a G(r/t), it is 
easy to see that C(r,t) = t a+1 C{r/t). Indeed, Fig. [3] 
shows a numerical evidence for the scaling of C, with 
a w 0. (If we require t — to be the origin, the scaling 
law is only accurate to 10 %. If we fit the location of the 
origin of time and drop the initial evolution (f < 200), 
this we get an accuracy of 3%. The displacement of the 
origin fits to about — 26 ). 



III. HARTREE APPROXIMATION OF SCALAR 
FIELDS 

In this section we review the inhomogeneous Hartree 
approximation and its application to our model. The 
reader can find a more detailed introduction in Ref. [3^ | . 
Its other name, Gaussian approximation, reflects the 
essence of the truncation of the dynamics: we disregard 
any connected higher n-point functions. Note in the con- 
text of the TV-component scalar field the leading order in 
1 /N expansion leads to very similar (also Gaussian) , but 
inequivalent approximation (44|. 

The operator equation in Heisenberg picture that we 



4 




0.5 1 1.5 2 

rescaled length (r/t) 

FIG. 3: Scaling of the equal time correlation function. Notice 
that the function has an approximate Gaussian shape, which 
is due to the nature of spinodal instability that created the 
domains. (There are deviations at small distance.) (L — 
1000, 1024 runs averaged) 



have to solve reads: 

(d 2 + m 2 )0(x,t) + ^ 3 (x,t) =0. (5) 
6 

We split off the quantum expectation value: (f>(x, t) — 
®{x,t) + tp(x, t), with (<p(x, t)) = 0. 

The Hartree approximation basically means that we 
restrict the density operator of the (p degrees of freedom 
to Gaussian at all times. For such degrees of freedom tp i 
we have the following identities: (<pf) = and wf(f>jj — 

We can simply take the quantum average of Eq. (|5j). 
or multiply from the left by 0(y, t y ) and get an equation 
for the Wightman propagator G < (x, t x ; y, t y ) = G < (x, y) 
by averaging that, too: 



dl + m 2 + ^G<(x,x) 



d 2 x +m 2 + ^ 2 (x) + ^G<(x,x) 



$(x) + ^ 3 (x) = 0,(6) 
G<(x,y) = 0,(7) 



It is remarkable that this simple truncation of the hier- 
archy of n-point functions leads to a self-consistent set of 
equations. Indeed, these are the Schwinger-Dyson equa- 
tions for the propagator in the leading order truncation 
of the two-particle irreducible (2PI) effective action [4|| . 

Gaussianity also implies that the Heisenberg operators 
at a finite time relate to the initial operators by a Bo- 
golyubov transformation: 



ip(x,t) 



d d k 
(2^p 



(8) 



So that we have the text-book operator at t — we 



set t/>e(a?,0) 



x 2wfe and ip^(x,0) = —iuk^t{x,Q). 



Here d is the number of space dimensions and uii = k 2 



with a contribution from the background field. The lad- 
der operators obey the usual [a- , at ] = (2n) d 5(ki — k-i) 
commutation relation. The initial particle spectrum are 
- ; — rA. These numbers appear in the 

k 

equal time two-point function: 



given by {aj.a% 



G < (x,t,x,t) = 



d d k 

(27T) 



1 \^(x,t)\ 2 (2nl+l) (9) 



This obviously diverges even in two space dimensions. 
The initial infinite mass shift we compensate by a mass 
renormalisation and introduce the finite mass squared m 2 
in the equation for ijj^: 



d d k 
(27r) d 



dl + m 2 r + -${x,t) 



[|^,t)r-^,0)| 2 ](2n| + l) (10) 



ipp(x,t) = 



A coupling renormalisation is also necessary in three 
dimensions [461 ] . but in this simple case we do not need 
to go beyond mass renormalisation. 

The traditional way of solving the dynamics in Gaus- 
sian approximation involves Eqs. ©, © and (TlTJ]) . One 
normally discretizes the equations on a space lattice. A 
consistent Bogolyubov transformation requires that the 
k index of the mode functions runs in the entire Fourier 
space of the lattice. In addition to the trivial background 
equation on an A 2 lattice this means A 4 complex equa- 
tions. 

But mode function expansion is just one of the possi- 
ble ways of solving Eqs. ([6]) and ([7]). Alternatively, we 
consider an ensemble of A e classical trajectories ipi(x,t), 
solutions of the equation 

dl + m 2 + ± [$ 2 (f, t) + (p 2 (x , t)) E ] ) <pi(2, t) = . 

Here (-} E stands for the ensemble average. Indeed, mul- 
tiplying the equation with ipi(y,t y ) and averaging over 
i (ensemble average), will bring us back to Eq. ((7|) with 
G K — > G e = (ip(x,t x )(f(y,t y )) E . But there will be no 
exact equivalence between ensemble and quantum aver- 
ages: the quantum two-point function G < is complex, 
G e is real. Notice, however, that the imaginary part of 
G < entirely decouples in Eq. ([7]), since the equal time 
propagator is always real. 

Of course, tfi(x, t) must be properly initialized to form 
a Gaussian ensemble of the correct standard deviation: 



m 2 , where m r 



is the renormalized mass, often amended 



(<p(x,0)<p(y,0)) 



d d k 
(2n) d 

d d k 
[2TT) d 



j-ik(x-y) 



(12) 
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Technically, we initialize tp in momentum space by a ran- 
dom phase and amplitude at the t — and t — St time 
slices. 

We intentionally introduced the factor h, as a control 
parameter for the fluctuation (p. This way we can tune 
strength of the back reaction of the quantum fluctuations 
to the background. In the classical theory it was possible 
to scale out A, here rescaling the field with 1/vA would 
also require to rescale Ti with A. If we stick to A = 6M 
in the numerics, it is the h in Eq. (|12[) that one can use 
to vary the coupling, effectively. 

Numerically, it is much simpler to solve N e N 2 real 
equations than A^ 4 complex ones, we found that even an 
ensemble of 1 <C N e <C N 2 was big enough. The simple 
structure of Eq. pip allows high speed implementations 
[55l ] . We note that the equation is not stable without 
manually fixing (ifi(x)) — after every leap-frog time 
step. 

Before embarking into the analysis of numerical re- 
sults, let us pause to discuss in what sense the Hartree 
equations represent a quantum correction to the classical 
dynamics. 

Notice that we can arrive at Eqs. ((6]) and (0 also from 
a different concept. Let us start a number of classi- 
cal trajectories from an initial Gaussian ensemble (e.g. 
Eq. (Tl2|) ). Instead of following the individual trajectories 
we can write down the equations for the n-point func- 
tions. Simply discarding the three or higher order cor- 
relators we get a closed set of equations, that coincide 
with Eqs. ([6]) and ([7]). We would also arrive to the same 
equations by truncating the 2PI effective action for the 
classical (or quantum) field theory to leading order. 

Indeed, whether we start from a classical or quantum 
Gaussian ensemble, the genuine quantum features start 
to appear if we keep the four-point equation at least. A 
self-consistent set of equations follows from the next-to- 
leading order truncation of the 2PI effective action, where 
one easily identifies the term, responsible for quantum 
effects jig. 

This statement, however, means that to Hartree or- 
der it is only the initial condition that reflects quantum 
physics. Do the mode function equations (|10|) or the 
propagator equation ([7]) introduce quantum corrections 
at all? 

The answer is yes. If we consider one single classical 
trajectory $, switching on h in Eqs. (TT2"]) will definitely 
enable many quantum phenomena, such as vacuum par- 
ticle production. Instead of doing Hartree, one can, of 
course, consider an ensemble of $ fields, initialized (as 
usual) with the just-the-half rule (analogous to Eqs. [T2"l) 
and evolve them classically. This classical ensemble will 
equally enable the same quantum phenomena, but it will 
bring in several classical artefacts, too, such as the decay 
of the quantum zero-point energy. These artefacts can 
most simply eliminated by shutting down all higher loop 
diagrams, down to the order where quantum and classical 
approximations agree: this is the Hartree approximation. 

Although it is possible to properly include higher or- 



der corrections [32j, |33[, they are not inevitable in the 
following two extremes: If the particle numbers are low, 
the higher order quantum corrections, that account for 
scattering of the quantum fluctuations, are not very im- 
portant compared to the dynamics of other energetic ob- 
jects, such as defects. If the particle numbers are very 
high, higher order quantum corrections are crucial, but 
they can be well estimated by a classical ensemble, here 
dominates in + 1/2 and the classical artefacts will 
then be suppressed. 

In the application considered in this paper we work 
with low particle numbers produced by the sparse net- 
work of defects. We believe that we do not miss the 
magnitude of the particle's back-reaction by ignoring the 
thermalisation of their spectral distribution. This semi- 
classical approach, however, approximates the damping 
of the classical degree of freedom at the lowest order. 



IV. CLASSICAL VERSUS HARTREE 
DYNAMICS 

The initial conditions given in section III Al define a 
(highly non-Gaussian) ensemble of domain wall config- 
urations at t = 15. For each member we define a 
(Gaussian) sub-ensemble of fluctuations. We follow the 
dynamics of this sub-ensemble in the Hartree approxima- 
tion. The final averaging over the domain wall configu- 
rations occurs at the very end of the calculation. At the 
time we switch on the quantum equations (|11[) we renor- 
malize the mass and thereby allow a smooth transition 
to quantum evolution. 

In Fig. 2] we show the evolution of the power spectrum 
of the background field. In classical field theory this is the 
only degree of freedom, whereas in the Hartree approxi- 
mation energy may drift into the "modes" (the ensemble 
of quantum fluctuations). 

The correlation length in Fig. [5] is defined by a Gaus- 
sian fit to the correlation function shown in Fig.[3J In har- 
mony with Fig. [5] we see no impact of the quantum fluc- 
tuations on the evolution of the macroscopic degrees of 
freedom. If all the domain wall loops were macroscopic, 
this would suggest that inverse total length of domain 
walls shown in Fig. [5] receives no significant quantum 
correction. Indeed, we again find a linear scaling, and 
we could not find a significant correction to the slope pa- 
rameter for t > 150. But the inverse domain wall density 
is not entirely linear. In the first half of the evolution it 
drifts away the classical solution. This reflects a transient 
decay of some classically more stable structures. 

The power spectrum in Fig.[4]does not give account on 
the created particles. The power spectrum of the ip x (x, t) 
functions in Eq. [TT] reflect the created particles. This 
spectrum does not scale, and performs a "boring" evolu- 
tion: only the amplitude changes slightly and always re- 
sembles the vacuum power spectrum. This confirms the 
assumption that the particles are created on the mass 
scale and not e.g. in the infrared. At and beyond the 
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FIG. 4: The scaling of the power spectrum. The macroscopic 
part (domain walls) clearly follows the scaling law both in the 
classical (top) and in the quantum (bottom) case. The scaling 
breaks at the tail of the spectrum (small classical structures) . 
(L — 500, average of 16000 and 192 runs for h = and 1, 
respectively.) 



mass scale the quantum fluctuations dominate over the 
classical structures in the power spectrum. 

By construction of the initial conditions the energy 
density transferred from the domain walls to the fluc- 
tuations does not significantly raise the temperature and 
so the thermal mass. It was an important assumption in 
our analysis that the wall width M is constant in time. 

The domain wall density and the correlation length are 
the key observables when we discuss scaling. Irrespec- 
tively to the coupling (A) or the strength of the fluctua- 
tions (h) we fit £ 3.4(2) - t. For the domain wall density 
we find C = l/L 2 as 1.66(3) /t for our initial condition, 
where I is the total counted length of domain walls on 
the lattice at a given time. It is remarkable that these 
dimensionless coefficients are robustly insensitive to the 
variation of the coupling or the value of ft. Also it does 
not depend on the lattice spacing nor on the initial noise 
amplitude or the details of the cooling procedure. 

To gain more insight into the small discrepancy be- 
tween the quantum and classical domain wall density we 
count the number of domains, and use this number to 
estimate the loop number density (n(t) = N(t)/L 2 ). We 
applied a cluster algorithm on the lattice and plotted 
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FIG. 5: The scaling of the correlation length (top) and the 
domain wall density (bottom) 



the resulting number density in Fig. [6l The effect of the 
quantum fluctuations is now striking. 

In the approximated quantum theory we get what we 
think we should: if the correlation length £ ~ t scales 
linearly, any number density must scale as n q ~ t~ 2 . 
The conclusion from Fig. [6] is that in the classical theory 
the domain number is dominated by microscopic struc- 
tures. One hint for the smallness of these "mini-domains" 
is that the inverse wall width M must appear in the 
n c i(t) ~ M/t scaling rule for dimensional reasons, and its 
coefficient is not extreme. We can directly measure an av- 
erage defect loop size by counting the loops (or domains, 
practically) that are in excess in the classical solution. 
The total wall length is also bigger in the classical case 
than with quantum correction. We used their quotient to 
estimate the size of these loops in Fig. [6j Unfortunately 
our numerics is not conclusive at later times, we expect 
this ratio to settle at a positive value ~ M _1 . 

To find out more about these small classical structures 
let us compare the lattice snapshots taken from the same 
run with and without quantum fluctuations. We picked 
the time t = 40 and cropped a larger lattice appropriately 
so that we can show the most phenomena in one image: 
Pig-El 

In these images the black regions correspond to the 
$ > domains. We plotted the 4> < domain with bet- 
ter colour resolution. The left hand side plot corresponds 
to the classical evolution, it shows many ripples and nu- 
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FIG. 6: Number of domains over lattice volume (top) in 
the classical and quantum framework. The classical scaling 
is counter intuitive. Dividing the difference of the total loop 
length by the number of loops (estimated by the number of 
domains) we get an average loop size (bottom). 



merous dark spots where $ goes locally close or beyond 
zero. In the middle (quantum) plot we see fewer ripples 
and some (but not all) of the dark spots of the classical 
plot are missing here. To the right, we demonstrate the 
excitation of the quantum modes by plotting (ip 2 (x, t)) E . 
These fluctuations are the strongest on the domain walls, 
which we interpret as particles in the bound state. 

Looking at a sequence of such snapshots gives more 
detail about how these ripples are produced. As the do- 
main walls shrink they emit classical waves, with a wave 
length of few times the wall width. These waves are 
equally present in the quantum case as well, where they 
are more damped. In the quantum fluctuation plot we 
also find traces of the classical ripples, but the spatial dis- 
tribution of the produced particles appears to be smooth, 
and the ripples in (i J 9 2 (a ? , t)) _ are an order of magnitude 
smaller than in $ 2 . 

On the snapshots we marked the interesting places by 
letters. In both sides of the letter "A" the ripples are lo- 
cally so high in amplitude that these spots are counted as 
a domain by the cluster algorithm and they contribute 
to the total length of domain walls. But they are not 
counted into the statistics of walls in the quantum case, 
since then their amplitude is within the threshold of zero 



field value. The amplitude of these spots actually oscil- 
lates, this is why we do not see the one on the right hand 
side of the letter "A" in the quantum plot. "B" marks the 
centre of ripples emitted earlier by the collapsing bubble 
marked with "D" . These are mostly damped in the quan- 
tum run. The waves around the bubble "D" are higher 
in amplitude than in the quantum case. Finally, there is 
a spot with strongly oscillating amplitude, marked with 
"C" . The magnitude of the quantum fluctuations oscil- 
late coherently with the background field value. 



V. DISCUSSION 

Let us summarize the numerical findings: The correla- 
tion length, which is fitted from the correlation function 
in direct space, reflects the macroscopic evolution. We 
find that the known scaling behaviour is unperturbed 
by quantum effects. On the microscopic level, however, 
where the scaling is broken, we find stronger quantum 
effects, as expected. 

We also find, that there are "mini-domains" in the 
classical simulation, that (at least partly) disappear if 
we switch on the quantum degrees of freedom. Its sim- 
plest explanation is that there are classically stable small 
structures that decay in a quantum field theory Now we 
can speculate what these could be. Natural candidates 
are oscillons, localised oscillating wave packets, which 
are (quasi-)stable solutions of the classical field theory 

If these small structures are indeed oscillons, their sta- 
bility is enhanced by low dimensionality. If in three di- 
mensions they are subject to a more rapid decay [50l.[5l|. 
making the quantum decay channel less significant and 
hence the quantum correction to the scaling even smaller. 

Indeed, a closer look on the lattice field revealed that 
there are small regions (with a diameter of 0(5) do- 
main wall width) that oscillate with a frequency ~ M. 
But oscillons are not the only structures that appear. 
The shrinking and collapsing domain walls emit classical 
waves with a wave length ~ M . We see these waves 
on the lattice snapshots as circular ripples. These ripples 
from various sources interfere and at the points of con- 
structive interference the field value may locally exceed 
zero and will then be counted as a small domain. 

Classical waves are emitted in the quantum field the- 
ory, too. In quantum mechanics this classical excitation 
is known as coherent state, which transforms into an en- 
larged width of the wave function, or into particles in field 
theory language. This is the point where quantum cor- 
rections enter: the classical waves are damped and their 
interference results in fewer and less stable localised os- 
cillating wave packets. 

In this picture there is a non-perturbative classical 
mechanism that converts the energy stored in the string 
(or domain wall) to microscopic objects. In a field theory, 
these objects are neither loop fragments, nor particles, 
but coherent oscillations of the field expectation value. 
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FIG. 7: Snapshot from a classical (left) and the corresponding quantum (middle) run. To the right the particle content is 
characterized by (^ifi 2 (x,t)^ . The darker points mean higher value. (These images were taken at t = 40 on a L = 128 lattice. 
We cropped out a piece of 75 x 75) 



Our numerics suggests that the scale of these classical 
waves are on the microscopic scale M. We observe that 
these waves are emitted from structures of size £, present- 
ing us with the challenge of explaining energy transport 
over a huge scale separation, Ml ~ 10 3 at the end the 
simulation. 



It is clear from the shown numerics that the domain 
wall decay was not enhanced by the quantum fluctuations 
and this conclusion we checked to stay true with fi = 2 
or A = 12. There is no indication for a direct decay chan- 
nel into particles. A direct decay might also manifest in 
the sensitivity to the choice of the lattice spacing as we 
switched between aM = 0.5 and 0.7, but we found no 
significant difference. However, the decay of the classi- 
cal waves and oscillons is no longer protected by scale 
separation. 

Finally, let us attempt to understand Fig. [5] The 
energy density associated to macroscopic _D-dimensional 
defects in d dimensions is ~ M 1+D t D ~ d . Their decay re- 
leases energy at a rate of ~ M 1+D t D ~ d ~ 1 . This energy is 
used to produce high amplitude classical structures (e.g. 
oscillons) that may have been counted as small domains. 
Since they emerge on the microscopic scale, their number 
density has a source of CsourceM ^^ 4 ^ 1 , where C source 
and the other constants we introduce here are dimension- 
less numbers of 0(1). These small structures can decay 
in various ways: a) In the quantum calculation we include 
the direct quantum mechanical decay into particles with 
a rate of T ~ M; b) The small objects can be hit by a 
domain wall or string, its rate is proportional to the de- 
fect density: C de f ect M D ~ d+1 t D ~ d ; c) These objects can 
also hit each other and annihilate. The probability of a 
given small object to meet an other one is proportional 
to its density n, which gives a rate of C co ;/M 1_d n. These 



together give the following equation for the density n 

M D-d+l M D 
fl+Tn+C 'defect — Jd^D — n + 

Ccoll-M Tl — Csource ^—D-\- 

. . (13) 

If the quantum decay into particles dominates giving a 
finite life time to these small classical structures, the den- 
sity n simply follows the source. Indeed we see n ~ t~ 2 
in Fig. [6] In the absence of T, however, we find that 
n ~ M/t solves Eq. (fT3| in consistence with our observa- 
tion. Since in this case n shows the same scaling as the 
domain wall density, counting them as defects does not 
spoil the observation of scaling. The classical approxi- 
mation of Eq. (|13p suggests that for d > 2 the collision 
term would dominate, giving n ~ t~ 2 . In higher dimen- 
sions, however, oscillons and other analogous structures 
are less stable, which introduces a decay term of classical 
nature bluring difference between classical and quantum 
scaling. 

VI. CONCLUSION 

In this paper we integrated the classical field equa- 
tions as well as the Hartree approximated quantum evo- 
lution of a scalar field in the broken phase, starting from 
a network of domain walls. The scaling of macroscopic 
observables was manifest also in the quantum theory. 
Our numerical results suggest that the direct decay of 
domain walls into particles is insignificant, as the per- 
turbative estimates predict. We can instead attribute 
the decay to the emergence of classical waves and other 
structures, such as oscillons. Since these coherent exci- 
tations of the quantum field theory are produced at the 
microscopic scale, their perturbative decay is no longer 
suppressed by the separation of scales. The production 
of large-amplitude classical oscillations is a genuine non- 
perturbative phenomenon that deserves further investi- 
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gation, as a similar effect is seen to drive the decay of 
cosmic strings in three dimensional field theory simula- 
tions 22]. Understanding the dominant decay channel 
of strings is of crucial importance for computing their 
observational signals. 
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